Topological qubits in graphenelike systems 
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The fermion-doubling problem can be an obstacle to getting half a qubit in two-dimensional 
fermionic tight-binding models in the form of Majorana zero modes bound to the core of super- 
conducting vortices. We argue that the number of such Majorana zero modes is determined by a 
Z 2 x Z 2 topological charge for a family of two-dimensional fermionic tight-binding models ranging 
from noncentrosymmetric materials to graphene. This charge depends on the dimension of the 
representation (i.e., the number of species of Dirac fermions - where the doubling problem enters) 
and the parity of the Chern number induced by breaking time-reversal symmetry. We show that in 
graphene there are as many as ten order parameters that can be used in groups of four to change 
the topological number from even to odd. 
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I. INTRODUCTION 

A major hurdle in the realization of quantum comput- 
ers is the problem of decoherence. Qubits generically do 
not last long in the presence of the environment. Over- 
coming decoherence is possible if the qubit is stored non- 
locally using a many-body state, in such a way that the 
reservoir, which only couples locally to the system, is un- 
able to damage the quantum information. tl An implemen- 
tation of this scheme can be achieved if the many-body 
ground state supports excitations obeying non-Abelian 
braiding statistics. Non-Abelian braiding statistics that 
departs from the Bose-Einstein or Fcrmi-Dirac statistics 
can only be realized in effectively two-dimensional (2D) 
systems, such as the m = 5/2 fractional quantum Hall 
state on the one hand,BI or in chiral p x ±ip y 2D supercon- 
ductors where a half vortex binds a zero-energy midgap 
state on the other handBa In the latter example, due to 
particle-hole symmetry (PHS), this zero mode is a Ma- 
jorana fermion. A two-state system - a qubit - can be 
assembled from one complex fermion made up of a pair 
of such Majorana fermions sitting at far away vortices. 
The splitting of the energies from exactly zero is expo- 
nentially small in the separation between the vortices and 
thus the time to degrade a quantum superposition is ex- 
ponentially large in the distance between the vortices. 

Majorana fermions bound to the core of vortices in a 
superconductor were discovered by Jackiw and Rossi.El 
There is a single zero-energy midgap state when the vor- 
ticity is n = ±1 (there are generically \n\ zero modes). 
However, these results are obtained for the case of the 
minimal representation of Dirac fermions whose support 
in space is the 2D continuum. The fermion-doubling 
problem, discovered in the context of lattice-gauge the- 
ory,Q prevents importing these results of Jackiw and Rossi 
to condensed-matter systems, which arc models defined 
on lattices. In graphene, for instance, one does have the 
Majorana fermions from which it is possible to assemble 
the qubits, as shown by Ghaemi and Wilczek.Q However, 
there arc four of them in each vortex because there are 



two Dirac cones for each spin polarization in graphene. 
Even numbers of Majorana fermions are not stable, as 
local perturbations can move them away from zero en- 
ergy. 

One can get much insight into the problem of how 
many Majorana fermions can be realized in effectively 2D 
tight-binding models if one looks into ideas for address- 
ing the fermion-doubling problem in lattice-gauge theo- 
ries. There is the original proposal due to Wilson, which 
is achieved by adding perturbations (Wilson masses) to a 
lattice Hamiltonian that opens gaps at undesirable dupli- 
cate Dirac points!! However, from the point of view of a 
lattice-gauge rcgularization of quantum chromodynam- 
ics, this option has the undesirable property of break- 
ing the chiral symmetry. Alternatively, the idea that 
the fermion doubling can be overcome by considering 
an n-dimensional system as a boundary of a (n + 1)- 
dimensional one was put forth by Callan and Harvey 
in Ref. ^| (see also Refs. [l(] and [ll]). In fact, it is the 
Callan-Harvey effect that is at wotk in the remarkable 
results obtained by Fu and Kane:ll3 (i) surface states 
in 3D topological insulators with time-reversal symme- 
try (TRS) realize an odd number of Dirac fermions in 
the minimal representation of the Clifford algebra, (ii) 
They can be used to achieve an odd number of Majorana 
fermions bound to vortex cores induced by the proximity 
to a typc-II superconductor. 

We argue in this paper that the Wilson prescription is 
a route to attain an odd number of Majorana fermions 
in effectively 2D condensed-matter systems. The chief 
reason is that one is not constrained to impose chiral 
symmetry as in lattice-gauge theory. However, not any 
Wilson mass can be used for this purpose, only those that 
break TRS. This approach naturally leads to a Z 2 x Z 2 
topological charge that discerns whether the system has 
an even or odd number of Majorana fermions attached 
to a superconducting (SC) vortex. In essence, the parity 
of the number of zero modes is determined by the num- 
ber of Dirac points which have not been knocked out by 
changing the Chern number via the TRS-breaking Wil- 
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FIG. 1: (Color online) Energy levels for the lattice Hamil- 
tonian (|s]) as a function of h- A = h z for h 1 = h A = 0.55h 3 , 
h 2 = — 0.45ft 3 , and a = A(r = oo) = 1. Although finite- 
size effects prevent the closing of the bulk gap, the vanishing 
curvature of the gap in the bulk is a footprint of each thermo- 
dynamic transition as a function of h z . In region IV, there are 
8 (=4x2) midgap states, four bound to the vortex and four 
bound to the antivortex. Their degeneracy is lifted by the 
combined effects of intravortex-level or intervortex-level repul- 
sion. There are 6 (=3x2), 4 (=2x2), and 2 (=1x2) midgap 
states in regions III, II, and I, respectively. In the limit in 
which the vortex and antivortex separation goes to infinity, 
the intervortex-level repulsion is exponentially suppressed and 
midgap states can be converted to Majorana fermions. Thus, 
an odd number of Majorana fermions attached to an isolated 
vortex are found in regions III and I. 



son mass. For systems where the number of species is 
odd, like in the case of surface states of TRS topologi- 
cal insulators, odd numbers of zero modes occur without 
breaking TRS (odd x even case). In the case of graphene, 
which is one focus of this paper and where there is an 
even number of Dirac cones, TRS must be broken so as 
to obtain an odd Chern number and, in turn, an odd 
number of Majorana fermions (evenxodd case). In all 
other cases, including graphene when TRS is unbroken 
(even x even case) and surface states of topological insu- 
lators with large enough magnetic field (odd x odd case), 
there are even numbers of Majorana fermions. Notice 
that, according to this Z 2 x Z 2 classification, systems 
defined on 2D lattices must have both SC pairing corre- 
lations and a nonzero Chern number that accounts for a 
thermal Hall effect in order to have non-Abelian quasi- 
particles. This is "the poor cousin," i.e., the mean- field 
version, of the v = 5/2 quantum Hall state. 



fermions in a 2D tight-binding model, we begin with the 
pure Rashba kinetic energy 



H„ := 



He 



(1) 



n— x,y 



that we define on a square lattice with sites denoted by 
r = rax + ny, where m and n are integers. Here, c\, = 
(cj. s ) is a doublet that creates on site r an electron with 
the spin projection s =j, J, along the quantization axis, 



cr 1 and <t~ = a 



<7 2 are the first two 
Pauli matrices while the third Pauli matrix er 3 defines 
the quantization axis in spin space, and the real-valued 
number a sets the energy scale for the Rashba hopping. 
At half filling, i.e., at vanishing chemical potential, the 
Fermi surface collapses to the four nonequivalcnt Fermi 
points 



p F = (0,0),(0 J 7r) ) (7r,0),( 7 r J 7r). 



(2) 



These Fermi points are TRS in that they change by a 
reciprocal wave vector under the inversion p — > p. Lin- 
earization of the energy spectrum of H a in the vicin- 
ity of these four Fermi points yields an 8 x 8 massless 
Dirac Hamiltonian, i.e., a reducible representation of the 
Clifford algebra four times larger than the minimal one 
in 2D continuum space. This is a manifestation of the 
fermion doubling. Hamiltonian (0) preserves TRS but 
breaks completely SU(2) spin-rotation symmetry (SRS). 

We now introduce the same spectral gap at all the 
Fermi points, Eq. (^) . We achieve this with the help of a 
singlet SC order parameter parametrized by a complex- 
valued A, 



H > 



E 



[a4(^) c 



H.c. 



that we add to Hamiltonian (1) 



(3) 



(4) 



This 2D tight-binding Hamiltonian is the relative to a 
noncentrosymmetric superconductor in the Rashba-Dirac 
limit and with singlet SC pairing studied in Ref. [l3|. A 
TRS-breaking vortex 



A(r)exp(i6) 



(5) 



with the profile A(r), where r and 9 are polar coordi- 
nates, binds four Majorana fermions at the Fermi energy 
via the Jackiw-Rossi solutions.E3 

Finally, we define the three independent TRS-breaking 
Wilson masses 



II. TUNING THE NUMBER OF MAJORANA 
FERMIONS 

To illustrate how it is possible to use a Wilson mass 
prescription to change at will the number of Majorana 



(6) 



parametrized by the triplet of energy scales, 

h = (h 1 , h 2 , h 3 ) ee (h A , h^, h z ) = (h x , h y , h z ). (7) 
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Each Wilson mass breaks TRS and breaks SU(2) SRS 
down to U(l). 

We are going to show numerically on the lattice and 
analytically in the continuum limit that the Wilson 
masses h can be used to change the number of Majo- 
rana fermions bound to the core of the vortex, Eq. (g), 
in 



H := H„ 



(8) 



one by one from four to zero. 

Figure [j] displays the energy spectrum of Hamilto- 
nian (@) obtained from numerical diagonalization on a 
square lattice made of 39x39 sites. Periodic boundary 
conditions are imposed in the presence of a vortex and 
an anti- vortex with winding numbers ±1, respectively, 
that are as far apart as possible. A continuum of energy 
cigenstates is visible as are bulk-gap-closing transitions 
as a function of h 3 = h z . Midgap states are also visible 
although they are not located at the Fermi energy be- 
cause of level repulsion for states attached to the same 
defect and because of the finite separation between the 
two defects. Starting from four midgap states per iso- 
lated vortex in the thermodynamic limit, increasing h 3 
decreases their number by one after each bulk-gap-closing 
transition. 

The same results follow analytically after lineariza- 
tion of the spectrum around the Fermi points, Eq. (^). 
Indeed, linearization of Hamiltonian (JsJ) yields, in the 
Bogoliubov-de-Gennes single-particle representation, the 
16 x 16-dimensional block-diagonal Hermitian matrix, 



H :=diag (H^H^H^Hi) 
with the 4x4 Hermitian blocks 



1~La — 



p_ 

6 A 



P 




s, ' 



V s s 



-p 



-p 

'Ij 



(9a) 



(9b) 



whereby the units h = v F = 1 have been chosen, the 
complex notation p = p 1 + ip 2 is used for the momenta 
whereby x denotes the complex conjugate of x, and 

5 1 = 5 2 = A, 8 3 = S 4 = A, (9c) 
Vi =V(a,o) = 2 {h 3 + h 1 +h 2 ), (9d) 
% =V(a,n) = -2(/i 3 + /i 1 -h 2 ), (9e) 

% = ^(tt.tt) = 2 i h S ~ h l ~ h 2> , (9 f ) 

V4 = V(n,o) = - 2 ( h 3 -K + h 2 ) ■ (9g) 

The mathematical form of any of the four matrices Hj is 
the same as that studied in Ref. provided the SC order 
parameter A and each r]j are identified, respectively, with 
the Kekule and Haldane masses in Ref.[T§ Thus, we can 
immediately borrow and tailor some of the results from 
Ref. [Fi] to the present case. 

If all the r) j 's are zero and the SC order parameter has 
a single vortex with unit winding number, there are four 



Majorana fermions bound to it. As the magnitudes of 
the 1 77 - 1 's increase, there will be a phase transition every 
time that |r/ 3 -| = |A(r = 00) |. Any such transition is 
characterized by a decrease in the number of Majorana 
fermions by one unit and a corresponding change in the 
value of the Chern number by ±1, depending on the sign 
of rjj. Therefore, by changing the Chern number of the 
system by ±1 each time, one can knock out the Majorana 
fermions one by one. 

Alternatively, one could start from the dominant Hal- 
dane masses limit defined by \rij\ > \A(r = 00) | with 
j = 1,2,3,4. In this limit, the system sustains the ther- 
mal integer quantum Hall effeG|t— (IQHE) and supports 
four chiral Majorana fermions. One can then change 
h so as to cross successive quantum phase transitions at 
which any one of the r) 1 , . . . , rj A equals in magnitude the 
spectral gap controlled by |A(r = 00) [. As before, each 
time we cross a phase transition, the Chern number and 
hence the number of Majorana fermions at a SC vortex 
core changes. 

We note that the presence of a nearest-ncighbor- 
hopping dispersion 



e(p) = -2t{cos Px + cosp ) 



(10) 



with f C a is equivalent to adding a constant chemical 
potential 



Mj = <P 3 ) 



(11) 



for each one of the four Fermi momenta, Eq. (g). The 
effect of this term is to shift the gap closing condition to 



j 



Al 2 



PI 



(12) 



i.e., our results can be generalized to systems with 
quadratic dispersions and naturally explain the results 
found in Ref. BtI 



III. TUNING THE NUMBER OF MAJORANA 
FERMIONS IN GRAPHENE 

We are now going to demonstrate that the very same 
control on the number of Majorana fermions achieved 
with Hamiltonian (||) is also possible in graphene. We 
recall that in graphene, electrons with spin s =t, i hop 
on a honeycomb lattice that is made of two triangular 
sublattices A and B. The conduction and valence bands 
touch at the two non-equivalent points K ± located at the 
opposite corners in the hexagonal first Brillouin zone (see 
Ref. [lj for a review). Finally, to account for the possi- 
bility of a SC instability, Nambu doublets are introduced 
with the index p and h to distinguish particles from their 
charge conjugate (holes). Hence, after linearization of 
the spectrum about the Fermi points K±, this leads to 
a single-particle kinetic energy represented by a 16x16- 
dimcnsional matrix 



a 1 p 1 + a 2 p 2 . 



(13) 
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TABLE I: The ten mass matrices with PHS that anticommute with ai and cti and commute with the singlet SC masses 
M ReSSC and M ImSSC . Each mass matrix can be assigned an order parameter for the underlying microscopic model. The latin 
subindex of the order parameter's name corresponds to the preferred quantization axis in SU(2) spin space. Each mass matrix 
either preserves or breaks TRS, SRS, and sublattice symmetry (SLS). Each mass matrix can be written as a tensor product 
X M1M2M3M4 = p Mi ® s M2 ® <r M3 ® t M4 , where /0 Mi , s M2 , <r M3 , and r M4 correspond to unit 2x2 and Pauli matrices that act on 
particle-hole, spin-1/2, valley, and sublattice indices, respectively. 



Order parameter 


TRS 


SRS 


SLS 


Here 


y 


IQHE 


False 


True 


False 


-^IQHE 


-^3003 


ReVBS, 
ReVBS y 
ReVBS z 


False 
False 
False 


False 
False 
False 


True 
True 
True 


^ReVBS^. 
-^ReVBS^ 
A^ReVBS^ 


^3110 
-^0210 
-^3310 


ImVBS,, 
ImVBS,, 
ImVBS z 


False 
False 
False 


False 
False 
False 


True 
True 
True 


-^ImVBS x 
jWlmVBS^ 
-^ImVBS^ 


-^0120 
^3220 
-^0320 


NCCly 

Neel z 


False 
False 
False 


False 
False 
False 


False 
False 
False 




-^3133 
-^0233 
-^3333 



Here, a 1 and a 2 are two 16 x 16-dimensional Dirac ma- 
trices. 

It was shown in Rcf. |lj that there exists 36 distinct 
order parameters such that any one, when added to % D , 
opens a spectral gap. These order parameters are iden- 
tified by seeking all 16x16 matrices from the Clifford al- 
gebra that anticommute with H D . One complex valued 
order parameter is that for singlet superconductivity. We 
shall denote the two corresponding 16x16 matrices from 
the Clifford algebra by M RoSSC and M ImSSC and define 
the perturbation 

H A := AiM ReSSC + A 2 M ImSSC (14) 

that opens the spectral gap 2| A| with the complex- valued 

A = A 1 +iA 2 (15) 

parametrized by the real-valued A 1 and A 2 when added 
to % D . Next, we seek all 16x16 matrices from the 
Clifford algebra that (i) anticommute with H D and (ii) 
commute with % A . In this way, we find all ten TRS- 
breaking order parameters listed in Table [| that alone 
would open a gap in the Dirac spectrum if not for 
their competition with the gap induced by singlet su- 
perconductivity. Within this set of ten matrices one 
can form groups of at most four matrices that are 
mutually commuting and therefore can be simultane- 
ously diagonalized. Here, we choose the four-tuplet 
{ReVBS^, , ImVBSy , Neel z , IQHE} for concreteness but 
the results hereafter apply to any other such four-tuplet 
of commuting mass matrices among the set of ten. Our 



main result regarding graphene is the fact that 
n = p ■ a 4 A 1 M ReSSC + A 2 Af ImSSC 

+ m 1 MrcVBS^ + m 2 -^ImVBS y 

+ m 3 M N6ch 4 i] M IQHE (16a) 
is unitarily similar to Eqs. (^)aj) and ( pb| ) with 



£1,2,3,4 = A, 






(16b) 


i]i = —m 1 4 m 2 


4m 3 


+ V; 


(16c) 


j] 2 = m 1 — m 2 4 


m 3 4 


'/■ 


(16d) 


T]3 = m 1 4 m 2 - 


m 3 4 


'/■ 


(16e) 


rji = —m 1 — m 2 


-m 3 


4 rj. 


(16f) 



The phase diagram in Fig. |2j follows. 

There is a total of four SC pair-potentials that open a 
uniform gap at K ± in graphcncO One is the singlet SC 
pair potential, which we have discussed so far and the 
remaining three are all triplet SC pair potentials. For 
each such triplet SC mass, as in the singlet SC, there are 
four competing orders that commute pairwise and can 
be used, in principle, to knock out one by one Majorana 
fermions bound to the cores of isolated vortices. 



IV. CONCLUSIONS 

In summary, we have identified a mechanism to over- 
come the fermion-doubling barrier that can prevent the 
attachment of an odd number of Majorana fermions to 
the core of SC vortices in graphenelike tight-binding mod- 
els. This mechanism relies on aZ 2 xZ 2 topological charge 
that measures the parity in the number of Majorana 



5 



|A| A 




FIG. 2: Schematic phase diagram of the competition between 
the singlet SC (A), thermal IQH (7/), magnetic bond (m lj2 ), 
and Neel (m 3 ) orders near the Dirac point in graphene. Here, 
m i,2,3 are fixed while \n\ and |A| vary. When m 12 .3 = 0, 
there are three phases (two large \q\ phases with 77 > and 
77 < and one large |A| phase, according to Ref. n4). If we 
choose m 12 ,3 in such a way that m 1 +m 2 + m 3 , m 1 +m 2 — m 3 , 
-mi + m 2 — m 3 , and —m 1 — m 2 + m 3 are all different, there 
are 15 phases as we change 77 and A. Shaded (nonshaded) 
regions represent a phase with the odd (even) Chern num- 
ber (proportional to the thermal Hall conductivity divided by 
temperature). In phases that are adiabatically connected to 
a nonsuperconducting state (the horizontal axis |A| = 0), one 
can switch off the pairing without closing the gap. In these 
phases, the number of Majorana fermions is thus even. 



fcrmions attached to an isolated vortex and the use of 
TRS-breaking order parameters that compete with each 
other and with the SC order parameter to knock out one 
by one the Majorana fermions. In this surgical way, an 
odd number of Majorana fermions can be made to bind 
the vortices in a singlet SC order parameter, whereas this 
could only be achieved for the more elusive triplet SC or- 
der parameter in Refs. ^, ||, and [l^. This mechanism 
applies to graphene with superconductivity induced by 
the proximity effect, provided a way can be found to also 
induce and select from the remarkably large variety of co- 
existing and competing order parameters that graphene 
supports those with odd Z 2 x Z 2 topological charge and 
thus odd number of Majorana fermions attached to iso- 
lated vortices. 
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